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Abstract 

Digital holographic microscopy is a fast three-dimensional (3D) imaging tool 
with many applications in soft matter physics. Recent studies have shown 
that digital holograms can be fit to electromagnetic scattering solutions to 
obtain the 3D positions of isolated colloidal spheres with nanometer preci- 
sion and millisecond temporal resolution. Here we describe the results of new 
techniques that extend the range of systems that can be imaged through fit- 
ting. We show that an exact multisphere superposition scattering solution 
can be used to fit holograms of colloidal clusters containing up to six spheres. 
We also introduce an approximate and computationally simpler solution, Mie 
superposition, that is valid for multiple spheres spaced several wavelengths 
or more from one another. We show that this method can be used to fit holo- 
grams of several spheres on an emulsion droplet, and we give a quantitative 
criterion for assessing its validity. 
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1. Introduction 

In-line digital holographic microscopy (DHM) has emerged as a powerful 
and useful 3-dimensional (3D) imaging technique for characterizing colloids 
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and soft materials [H El El IH El El El • Unlike conventional microscopes, holo- 
graphic microscopes use a coherent source to illuminate a sample, generating 
a two-dimensional interference pattern that encodes 3D information (Fig- 
ure [T^). For relatively simple objects, such as colloidal spheres, the most 
precise way to recover this 3D information is to fit the hologram to an elec- 
tromagnetic scattering solution such as the Lorenz-Mie solution. The fitting 
technique can reveal the positions of particles with nanometer-scale spatial 
precision [1] on sub-millisecond time scales. The high temporal resolution, 
more than an order of magnitude larger than that of other 3D imaging tech- 
niques such as confocal microscopy, has enabled observations of new physical 
phenomena, as demonstrated in a recent study [7] showing the detailed dy- 
namics of a colloidal particle interacting with an oil-water interface. 

Thus far, the fitting technique has only been used to analyze holograms 
of simple scatterers. The seminal work of Lee et al. [1] as well as recent ex- 
periments from Kaz et al. |7] examined single isolated spheres, for which the 
well-known Lorenz-Mie solution suffices. Recently Fung et al. used a mul- 
tisphere superposition code [8] to fit digital holograms of non-rigid colloidal 
sphere dimers and rigid trimers [6]. To date, however, there have been no 
reports on using scattering solutions to analyze holograms of more complex 
scatterers, including clusters of more than 3 particles or colloidal systems 
with configurations that are unknown a priori. 

Such systems are of significant physical interest. Recently, Meng et al. 
showed how clusters of six or more spheres with a short-ranged attraction self- 
assemble into morphologies governed by rotational and vibrational entropy 
[9]. Detailed 3D observations of the self-assembly process and transitions 
between different cluster morphologies were not possible in that work, which 
relied on conventional microscopy. Another interesting system consists of 
colloidal particles bound to the surface of an emulsion droplet. The particle 
configurations, which need not be rigid or symmetrical pUj, may reveal infor- 
mation about the still poorly-understood interaction forces between colloidal 
particles on a curved liquid-liquid interface [TT] . 

In this work we demonstrate new methods to fit digital holograms of 
more complex scatterers, including rigid clusters containing up to 6 particles 
and nonrigid configurations of particles bound to the surface of an emulsion 
droplet. We used two analytical scattering solutions to fit our experimental 
data: an exact multisphere superposition solution and an approximate, com- 
putationally simpler approximation, Mie superposition. Both solutions can 
be used to compute holograms from an arbitrary number of spherical parti- 
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cles in an arbitrary geometry, but they differ in their precision to which they 
can fit holograms of the two experimental systems. We discuss the applica- 
bility of the Mie superposition approximation and the challenges remaining 
in the application of these techniques to other colloidal systems. 



2. Methods 

2.1. Recording Digital Holograms 

We recorded digital holograms on a modified Nikon inverted microscope 
(Figure[l|D). Our apparatus is described in detail in [51 [7]. In brief, we coupled 
a 660 nm laser diode (Opnext HL6545MG) into a single-mode fiber, from 
which we obtained up to 60 mW of laser power. The beam was collimated 
by the microscope's condenser before striking the sample. To image colloidal 
clusters, we used a 60x, 1.20 NA water immersion objective (Nikon). To 
image particles on droplets, we used a 100 x, 1.40 NA oil immersion objective 
(Nikon) with an immersion oil of index 1.412 (Cargille). In both cases, the 
index of refraction of the immersion medium was chosen to match that of the 
continuous phase of the sample, thus minimizing height-dependent spherical 
aberration [12j. Images were recorded by a Photon Focus MVD-1024E-160 
camera and captured to disk by an EPIX PIXCI E4 frame grabber. 

We also integrated an optical trap into the apparatus to manually as- 
semble colloidal clusters and manipulate them. The trap beam was kept off 
during holographic imaging. Light from the 830 nm trap diode (Sanyo DL- 
8142-201) was filtered through a single-mode fiber and reflected by a dichroic 
mirror into the back aperture of the Nikon objective lens used for imaging, 
as shown in Figure [T]d. 

2.2. Sample Preparation 

We imaged two different types of samples in this work: colloidal clus- 
ters bound by depletion forces, and particles bound to the surface of index- 
matched emulsion droplets. 

2.2.1. Colloidal Clusters 

We manually assembled tetrahedra containing four spheres, trigonal bipyra- 
mids containing five spheres, and polytetrahedra containing six spheres using 



the optical trap described in Section 2.1 Bright-field micrographs of these 



clusters are shown in Figure [2] To make the clusters, we first prepared a 
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Figure 1: a) Schematic illustration of in-line digital holographic microscopy. An incident 
plane wave illuminates a scatterer, and a detector records the interference pattern formed 
by the scattered light and unscattered incident light. The inset shows a typical interference 
pattern, or hologram, formed by a single 1-^m-diameter polystyrene sphere in water, b) 
Schematic illustration of DHM optical train, showing the counterpropagating laser beams 
used for DHM (lighter red) and optical trapping (darker red). 




Figure 2: Bright field micrographs of clusters made from 1.3-fim diameter polystyrene 
spheres, recorded with a 60x 1.20 NA water immersion objective. Ball-and-stick models, 
indicating cluster geometry and orientation, are a guide to the eye. Scale bar 5 fim. Ug 
denotes the number of spheres in each cluster, a) Tetrahedron (jig =4). b) Trigonal 
bipyramid (n^ = 5). c) Polytetrahedron (n^ = 6). 
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dilute suspension of monodisperse, 1.3-/im-diameter surfactant-free, sulfate- 
stabilized polystyrene spheres (Invitrogen) at a volume fraction of 8 x 10~^ in 
an aqueous solution containing 5 niM NaCl and 246 mM sodium dodecyl sul- 
fate (SDS). We loaded this suspension into sample cells made from untreated 
glass, as described in [B]. Because the SDS concentration was far above the 
critical micelle concentration [13], the SDS formed micelles which induced a 
short-ranged depletion attraction between the polystyrene spheres [m [15] . 
This attraction caused spheres to bind together into clusters when drawn into 
the focus of the optical trap. In all cases, we verified the cluster geometry by 
direct observation with bright-field microscopy before holographic imaging. 

2.2.2. Particles on Emulsion Droplets 

We prepared emulsion droplets laden with colloidal spheres by dispers- 
ing a suspension of microspheres in oil into an aqueous continuous phase. 
We suspended 0.8 /xm diameter polymethyl methacrylate (PMMA) spheres 
coated with poly(12-hydroxystearic acid) stabilizer [161 [H] in decane at a 
mass fraction of 2.1 x 10~^. The aqueous continuous phase contained 0.2% 
w/w Pluronic P123 triblock copolymer surfactant and 56% w/w glycerol. 
The glycerol matched the refractive index of the continuous phase to that of 
decane, so that only the PMMA spheres scattered light. We prepared the 
emulsions by mixing 0.5 mL of the PMMA-containing decane with 20 mL of 
the continuous phase in a 40 mL scintillation vial and shearing the mixture 
for 3 minutes at 9500 rpm with an Ika T9 Basic homogenizer equipped with 
a S25N-8G dispersing tool. After emulsification, the emulsion was diluted to 
17% v/v with additional continuous phase. A micrograph of a typical emul- 
sion droplet, laden with 4 particles but with the continuous phase slightly 
mismatched to allow the droplet to be seen, is shown in Figure |3j 

In order to minimize both the effects of nearby glass surfaces on the inter- 
particle interactions and unwanted back-reflections of the scattered light dur- 
ing imaging, we used negative dielectrophoresis (NDEP) to keep the particle- 
laden decane droplets away from the top of our sample cells, to which they 
would otherwise have risen due to buoyancy. We prepared sample cells sim- 
ilar in design to those used for the cluster measurements, except that the 
surface of the uppermost glass microscope slide in contact with the emulsion 
contained a linear array of interdigitated indium tin oxide (ITO) electrodes 
with a 40 /xm spacing between adjacent electrodes [HI [191 [20]. By applying 
a 10 V peak-to-peak, 300 kHz square wave with an arbitrary waveform gen- 
erator (Agilent AFG 3022B), we produced a downward NDEP force on the 
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droplets. Diffraction from the edges of tlie ITO electrodes was negligible. 

We used photolithography and wet etching to prepare the interdigitated 
electrode arrays on microscope slides coated with a 30 nm layer of ITO 
(Delta Technologies, CB-901N coating). We spin-coated the ITO surface of 
the slides with Shipley S1813 positive photoresist at 5000 rpm and soft baked 
the slides on a 115°C hot plate for 1 minute. We first defined the electrode 
pattern by exposing the photoresist to UV light through a photomask (150 
mJ/cm^ exposure at 405 nm), then developed the photoresist by immersion 
in Microposit MF CD-26 developer for 1 minute at room temperature. Fol- 
lowing an overnight hard bake in a 90°C oven, we etched away the excess 
ITO with an aqueous solution containing 40% v/v HCl and 10% v/v HNO3 
for 12 minutes at room temperature. 

2.3. Fitting Digital Holograms 

After recording digital holograms and suitable backgrounds in an exper- 
iment, we analyzed the holograms by fitting them to scattering models. We 
have released our software tools in the open-source package Holopy, avail- 
able at 'https : //launchpad.net/holopy, As described in detail in [6j , we 
divided the recorded holograms by a background taken in an empty field of 
view and normalized them to a mean value of 1 . We then used the Levenberg- 
Marquardt algorithm to minimize the sum of residuals between the exper- 
imental holograms and model holograms computed from a set of fixed and 
variable physical parameters, including particle coordinates, sizes, and re- 
fractive indices. We discuss the models in more detail in Section [H 

When fitting holograms of rigid clusters, we assumed that all of the 
spheres in the clusters had the same radius and real refractive index, and 
that the gap between spheres was uniform. We varied the common radius, 
the common particle refractive index, the common gap distance, the cluster 
center of mass, the three orientational Euler angles, and the scaling coeffi- 
cient a (see equation [T]) . When fitting holograms of particles on an emulsion 
droplet, we varied the common particle radius, the common particle refrac- 
tive index, 3 spatial coordinates per particle, and the scaling coefficient a. 

3. Modeling Holograms 

In our model, the normalized hologram Ihoi may be described as 

Ihol = l + 2a^ [Escat ■ e] + a^\Escat\^ (1) 
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where ^scat is the scattered electric field, and e is the unit vector describing 
the incident polarization. As described in [U |6] , a is inversely proportional 
to the magnitude of the incident electric field Ej„c and allows us to model 
variations in the incident laser intensity. In our experiments, we know the 
optical parameters of the system, including the medium refractive index, the 
incident laser wavelength, and the apparent detector pixel size (after mag- 
nification). Therefore, to model a hologram of some given configuration of 
particles, we need only an algorithm for computing Egcat- Any such algo- 
rithm must include the full radial dependence of Egcat on spherical Hankel 
functions of the first kind and their derivatives, as the far-field asymptotic 
approximation does not hold for the distances at which we record holograms 

We used two analytical methods to calculate Egcat for multiple spheres. 
The first method, described in detail in [6| , uses the multisphere superposition 
code SCSMFO developed by Mackowski et al. |8]. Previously we termed this 
approach the "T-matrix method" , but we now prefer the term "multisphere 
superposition" since we never actually compute the T-matrix itself, which is 
typically used for orientational averaging. This method accounts for multiple 
scattering and near- field optical interactions between the particles [HI 121] . 
The second approach is to superpose the fields calculated from the Lorenz- 
Mie solution for each of the spheres, taking into account the phase differences 
arising from the displacement of the spheres along the optical axis. This 
approach, which we term Mie superposition, assumes that only a plane wave 
falls onto each sphere. It neglects multiple scattering and near-field coupling. 
Escat calculated from Mie superposition is the lowest-order approximation to 
multisphere superposition; it is the initial solution that SCSMFO iteratively 
corrects in the order-of-scattering approach |22j . 

The main advantage of Mie superposition over multisphere superposition 
is that it is far less computationally costly for multiple spheres that are spaced 
far apart. SCSMFO expresses superposition of vector spherical 

harmonics whose origin lies at the center of mass of all the particles. It 
does so using vector spherical harmonic translation theorems to translate 
field expansions centered about each sphere to the cluster center of mass 
[Sj. When the particles are far from their center of mass, the final field 
expansion must contain many more terms than would be necessary in a 
field expansion centered about the individual spheres. We therefore use the 
Mie superposition approach for our submicron-diameter PMMA particles on 
emulsion droplets, which are weakly scattering (size parameter x ~ 5 and 
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relative index m ^ 1.07) and spaced far apart. We discuss the validity of 
Mie superposition further in Section |4.3[ 

We assess the goodness-of-fit between our recorded holograms and the 
best-fit models using two statistical measures. The first measure is chi- 
squared per pixel, x% the quantity the fitting algorithm attempts to mini- 
mize: 



TV 

2 _ ' " 



1 
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The sums run over all pixels of the recorded normalized hologram Ihoio 
and the best-fit model hologram If a. For any given hologram, comparing 
to an expected noise level allows us to assess whether deviations between the 
recorded hologram and best-fit model are due to instrumental noise or to a 
systematic error in the model's description of the underlying data. Assuming 
noise in the least significant bit of an 8-bit camera, we would expect values 
greater than (1/255)^ = 1.54 x 10^^ to originate from systematic errors. 
The second statistical measure we use is the coefficient of determination 
We define as 



] (Iholo Iholo) ] (Iholo 



where Ihoio is the mean value of the recorded hologram, which is 1 by our 
normalization [23] . measures the fraction of the variation of the recorded 
hologram from its mean value that is captured by the best-fit model, in- 
dependent of the amount of variation in the hologram. Unlike x^, R^ can 
be compared across physical systems that differ significantly in scattering 
cross section and hence hologram fringe amplitude. Thus R^ is more useful 
for assessing the goodness-of-fit of a best-fit model hologram when system- 
atic differences between the best-fit model and the recorded data dominate 
instrumental noise. 



4. Results and Discussion 

4.1. Rigid Clusters 

Figure |4] compares recorded and best-fit holograms, calculated with mul- 
tisphere superposition, for tetrahedral, trigonal bipyramidal, and polytetra- 
hedral clusters. Qualitatively, we observe excellent agreement between the 
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Water, glycerol, Pluronic P123 




Figure 3: Bright field micrograph of decane droplet with 4 PMMA spheres on its surface, 
recorded under 100 x magnification and differential interference contrast (DIG). Here, a 
slight index mismatch between the decane and the continuous phase makes the droplet 
visible under DIG. Scale bar 5 /zm. 



Cluster Radius (nm) Xp ^ 

Tetrahedron 670 ± 30 1.18 x 10^^ 0.923 

Trigonal bipyramid 640 ± 30 9.48 x 10"^ 0.910 

Polytetrahedron 650 ± 20 1.24 x lO'^ 0.877 

Table 1: Fitted radii and goodness-of-fit statistics Xp and for rigid clusters holograms 
in Figure |4] 
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Figure 4: Recorded and best-fit model holograms of rigid clusters. The third column 
compares the recorded holograms (solid lines) to the best-fit models (open symbols) along 
the color-coded dashed horizontal lines shown in the holograms. The renderings in the 
rightmost column show the cluster orientations determined from the fits. In the render- 
ings, the incident light propagates into the page, a) Tetrahedron {us =4). b) Trigonal 
bipyramid (n^ = 5). c) Polytetrahedron (n^ = 6). 
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recorded and best-fit holograms. In particular, the best-fit models reproduce 
the highly non-axisymmetric fringes in the recorded holograms, which de- 
pend strongly on the cluster orientations. The quality of the agreement is 
confirmed by the values of the fits, which are close to 1 (Table [T]). Also, 
the fitted particle radii are close to the manufacturer's reported value of 650 
nm. However, the values of Xp we observe are an order of magnitude larger 
than what we would expect due to camera noise, and indicate that further 
improvements to fits will depend on modeling additional physical phenomena 
or improving the convergence of the fitter. 

4-2. Particles on Droplets 

Here we show that our fitting techniques may also be applied to multi- 
sphere systems without a fixed geometry. Figure |5] shows a hologram of six 
particles bound to the surface of a decane droplet and a best fit model cal- 
culated using Mie superposition. Again, the qualitative agreement between 
the fringes of the recorded hologram and the best fit model is good. Quanti- 
tatively, we find Xp = 7.29 x 10~^ and = 0.811. Xp is much lower than the 
values we obtained for the clusters, primarily because the peak amplitude of 
the hologram in Figure |5] is significantly smaller than the peak amplitudes 
of the cluster holograms in Figure |4| The value of R^ indicates that the fit 
is slightly worse that the fits for the clusters. However, because we know 
the particles are bound to the surface of a spherical droplet, we can inde- 
pendently test the accuracy of the fitted particle positions. While we cannot 
directly image the decane droplet, which is index- matched to the continu- 
ous phase, we can fit the particle coordinates to the surface of a sphere, as 
shown in Figure |5]i. The average difference between the radial distance of 
each particle from the droplet center and the fitted droplet radius is 60 ± 60 
nm. Differences of this scale are comparable to previously reported preci- 
sions for DHM m ^ , and may be partially accounted for by variations in the 
interfacial contact angle between different particles [7]. 

To determine whether the slightly worse value of R^ obtained in this fit 
is due to the Mie superposition approximation, we examine the validity of 
this approximation in further detail. 

4-3. Applicability of Mie Superposition 

To confirm that Mie superposition is a suitable means for analyzing holo- 
grams like that in Figure [5} in which multiple weakly scattering spheres are 
situated several diameters apart, we fit the same hologram using multisphere 
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a) Recorded hologram b) Best fit c) Comparison d) Rendering 




50 100 150 200 250 
Hoiogram coiumn number 



Figure 5: Holograms of six PMMA spheres on a decane droplet, a) Recorded hologram, 
b) Best-fit model, c) Comparison between recorded hologram (solid lines) and best-fit 
model computed using Mie superposition (open symbols) along the dashed horizontal 
lines in the holograms, d) Rendering showing the sphere positions determined by fitting 
the holograms. The incident light propagates into the page, and the scale bar is 1 /im. 
The larger blue sphere is a guide to the eye; its position and diameter, 4.35 /im, were 
determined by fitting the coordinates of the six particles to a sphere. The small red 
sphere indicates the particle showing the largest discrepancy in position along the optical 
axis between fits to Mie superposition and multisphere superposition. 
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superposition. The multisphere superposition fit yielded Xp = 8-08 ^ 10~^ 
and = 0.790, comparable to the Mie superposition fits. Differences be- 
tween fitted particle coordinates in the two in-plane directions, perpendicular 
to the optical axis, were at most 22 nm. The largest difference between fit- 
ted coordinates along the optical axis was 154 nm. The sphere showing the 
largest displacement along the optical axis is highlighted in red in Figure 
|5ji. The size of these differences, as well as the lack of improvement in the 
goodness-of-fit using the multisphere superposition code, indicate that Mie 
superposition is an appropriate approximation. We conclude that the smaller 

values for this system stem from physical effects — such as weak scatter- 
ing by the decane droplet — that neither multisphere superposition nor Mie 
superposition account for. 

Further insight into the applicability of Mie superposition comes from 
examining the particle showing the largest coordinate difference along the 
optical axis. As shown in Figure [5]i, the largest difference occurs when two 
particles nearly occlude one another. When such a configuration occurs, the 
assumption that the field incident on each sphere is simply the illuminating 
plane wave is clearly invalid, as the colloidal spheres scatter most strongly 
in the forward direction. The field incident on the occluded sphere should 
therefore include a significant component of the scattered field from the first 
sphere. Whereas the multisphere superposition solution accounts for this 
multiple scattering effect, Mie superposition does not. 

As expected, Mie superposition fails to fit holograms from the clusters 



discussed in Section |4.1[ where the constituent spheres have a relative index 
m ^ 1.2 and size parameter x ~ 8.3. A fit to the polytetrahedron hologram of 
Figure Ep using a Mie superposition model yielded xl = 5.39 x 10^^ and R^ = 
0.463 (Figure [6]). These values are much poorer than the values associated 
with the multisphere superposition model. Moreover, qualitative differences 
between the best-fit Mie superposition model and the experimental hologram 
are readily apparent. 

To quantify the validity of Mie superposition for calculating holograms — 
or any other scattering-related physical quantity — of multiple spheres, we 
propose a dimensionless figure of merit. Mie superposition assumes that the 
exciting field at any sphere, Eg^, is approximately equal to the incident plane 
wave Ej„c- Following the multisphere superposition approach, consider Eex 
for any given sphere i to be equal to the sum of the incident plane wave and 
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the scattered waves from every other sphere at i: 

N 

Eex,j Ej^c -|- ^ ^ EgQQ^f j. (4) 

For Mie superposition to be vahd, |Escat,j| must be much smaller than |Ej„c|- 
This requires the particles to be far enough apart that their near fields do 
not couple. Then, |Escaij| scales approximately as \Eiinc\S/kR, where i? is a 
typical interparticle distance and S denotes the magnitude of the amplitude 
scattering matrix of sphere j in the Mie solution. From the optical theorem, 
in the Mie solution the extinction efficiency Q^xt for a single sphere is given 
by 

Qext = ^^[SiQ)] (5) 

where 5(0) is the complex amplitude scattering matrix evaluated in the for- 
ward direction j;24j. Then, if Qext^^/kR ^ 1, lE^catjl ^ |Ej„c| and Mie 
superposition should be valid. 

Taking the fitted drop diameter of 4.35 /im as a typical interparticle 
spacing for PMMA spheres like those in the droplet experiments, we find 
Qextx"^ /kR = 0.11, in agreement with our previous conclusion that the Mie 
superposition approximation is accurate for this sample. However, if we 
consider the sphere shown in red in Figure [5]i, using the nearest-neighbor 
distance of 1.19 fim for R yields QextX^/kR = 0.35. This larger value indi- 
cates that Mie superposition is a poorer approximation for this particle, as 
borne out by the 150 nm difference in the fitted position from Mie superposi- 
tion and multisphere superposition. In contrast, for the polystyrene spheres 
used in the cluster experiments, which are separated from each other by 
approximately a particle diameter, QextX^/kR = 13. We therefore conclude 
that QextX^ /kR should be approximately 0.1 or smaller for Mie superposition 
results to be trusted to a precision of 100 nm or better. 

When QextX^/kR <^ 1 and kR is large, Mie superposition is not only 
valid but advantageous compared to multisphere superposition. On a desk- 
top equipped with a 3 GHz Intel Core Duo processor, computing a 256 x 256 
hologram for the configuration of particles illustrated in Figure |5]i required 
11.8 seconds for multisphere superposition and 3.0 seconds for Mie super- 
position. In contrast, computing the 200 x 200 hologram of the six-particle 
polytetrahedron of Figure [4}:; required 3.9 seconds for multisphere superposi- 
tion and 1.8 seconds for Mie superposition. Here multisphere superposition is 
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only 2 times slower, an acceptable penalty given that the Mie superposition 
fit fails. 

4.4- Improving the Fits 

Our values of Xp are larger than the expected noise level. This indicates 
that our modeling and fitting procedures, rather than instrumental noise, 
limit the goodness of our fits. Challenges for future work fall into two cat- 
egories: modeling additional physical phenomena, and properly optimizing 
model parameters to determine a best fit. 

Our current hologram models do not capture several physical effects. 
First, while we choose immersion fluids for our microscope objective lenses to 
index-match the continuous phases of our samples as closely as possible, we 
cannot completely eliminate spherical aberration. Second, although we have 
fit for only a single, average particle radius in these experiments, no colloidal 
suspension is perfectly monodisperse. Third, while the well depth of the de- 
pletion attraction in the clusters is at least several fc^T, thermally-induced 
vibrations could perturb the cluster structures. We are currently working to 
model and account for spherical aberration in our hologram analysis. Fitting 
for individual particle positions in a cluster, rather than assuming geometric 
regularity, may address the latter two phenomena. 

Local minima in the fitting landscape pose another challenge to obtaining 
good fits. The number of local minima may increase with the number of 
parameters in the hologram models. This effect will become more significant 
when we study more complex scatterers or incorporate additional physical 
parameters, such as particle polydispersity, into our hologram models for 
multiple spheres. In order to prevent fits from becoming trapped in local 
minima, we are currently exploring the use of other minimizers, including 
one based on the r-algorithm [25] . 

5. Conclusions 

We have demonstrated that superposition solutions are accurate models 
for fitting digital holograms of up to six colloidal spheres. The multisphere 
superposition solution can reproduce the holograms of clusters in which there 
is significant coupling of the near-fields between spheres. We have also intro- 
duced a faster computational tool, Mie superposition, that should facilitate 
the study of multisphere systems where the particles are far apart. This 
approximate technique is accurate when QextX^/kR ^ 1. 
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The six-particle clusters we have imaged are large enough to show non- 
trivial dynamics related to multiple ground states [9]. We have also imaged 
configurations of spheres without a fixed geometry. Our experiments are thus 
the first step towards measuring internal rearrangements and self-assembly of 
sphere clusters. Moreover, analysis of dynamical information obtained from 
time series of holograms like those studied in Figure [5] may yield insight into 
colloidal interactions at liquid-liquid interfaces. 
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a) Recorded hologram b) Best fit, Mie superposition c) Comparison d) Rendering, 
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Figure 6: Fit of polytetrahedron hologram in Figure |4j: performed using Mie superpo- 
sition, a) Recorded hologram (shown again for ease of comparison), b) Best-fit model 
determined from Mie superposition, c) Comparison between recorded hologram (solid 
lines) and Mie superposition model (open symbols) along the dashed lines in the holo- 
grams, d) Rendering showing cluster orientation determined by Mie superposition fit. 
The incident light propagates into the page. 



19 



